Anisotropy of the Upper Critical Field and Critical Current in Single Crystal MgB 2 
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We report on specific heat, high magnetic field transport and ac— susceptibility measurments on 
magnesium diboride single crystals. The upper critical field H C 2 for magnetic fields perpendicular 
and parallel to the Mg and B planes is presented for the first time in the entire temperature range. 
A very different temperature dependence has been observed in the two directions which yields to a 
temperature dependent anisotropy with V ~ 5 at low temperatures and about 2 near T c . A peak 
effect is observed in the susceptibility measurments for fioH ~2 T parallel to the c— axis and the 
critical current density presents a sharp maximum for H parallel to the ab— plane. 

PACS numbers: 74.25.Bt, 74.25.Dw, 74.60. Ec 



I. INTRODUCTION 



Since the discovery of superconductivity in magnesium 
diboride at 39 K in January 200ltl enormous amount of 
work has been done which helped to elucidate many of its 
physical properties. Among others thcLjpultiband elec- 
tronic structure proposed theoreticallycra with the mul- 
tiple superconducting energy gaps has been approved. 
The two main superconducting energy gaps have already 



been experimentally evidenced by different techniques 
like for instance specific heat measurmentsQ or Andreev 
reflectionO. A larger gap is attributed to two-dimensional 
Px-y orbitals and a smaller gap to three-dimensional p z 
bonding and antibonding orbitals. Such a picture in- 
dicates a significant anisotropy of the superconducting 
state. Data reported so far on the anisotropy factor 
L = H C 2\\abl Hc2±ab scatter from 1.1 to 13 depending on 
the form of material (polycrystals, thin films, single crys- 
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tals) and method of evaluation. The problem remains to 
be clarified on high quality single crystals. 
The paper is organized as follows. The experimental 
details are given in section II. In section III, we show 
a consistent way of extracting the upper critical field 
H C 2 for magnetic fields parallel and perpendicular to 
the basal ab planes. H C 2 has been determined in the 
entire temperature range by high field magnetotrans- 
port, and those values are compared to the ones deduced 
from ac-susceptibility below 5 T and specific heat mea- 
surments below 7 T. The perpendicular critical field 
reveals a conventional character with a linear tempera- 
ture dependence near T c and HoH C 2± a b(0) — 3.5 T while 
H C 2\\ a b shows a positive curvature at temperatures above 
20 K with fJ.oH C 2\\ a b(0) ~ 17 T. As a consequence the 
anisotropy factor T is temperature dependent with T(0 
K) ~ 5 and T(near T c ) about 2. The angular depen- 
dence of the upper critical field measured at 5.4 K and 
26 K show an elliptic form as predicted by a one-band 
Ginzburg-Landau theory, but obviously this theory can 
not account for the temperature dependent T parameter. 
Finally, the magnetic field and angular dependence of 
the critical current deduced from ac— susceptibility mea- 
surments is discussed in section IV. 



II. SAMPLE PREPARATION AND 
EXPERIMENT 



Experiments have been performed on high quality 
MgB2 single crystals coming from one batcbB. The crys- 
tals show clear hexagonal facets with flat and shiny 
surfaces. The typical dimensions are 50x50x10 /im 3 . 
The resistivity has been measured at 13 Hz as a function 
of magnetic field up to 28 T for different temperatures 
and angles (6) between H and the a6-planes of the crys- 
tal. The magnetic field was always orthogonal to the 
measuring current. Gold electrodes have been evapo- 
rated as stripes overlapping the top plane of the sample 
with contact resistance of about 1 in the four-probe 
configuration. The ac— susceptibility was deduced from 
the local transmittivity measured with a miniature Hall 
probe. The ac excitation field (h ac ~ 3G, to ~ 23Hz) was 
perpendicular to the ab— planes and to the Hall probe 
plane and superimposed to a dc field which was making 
an angle < 9 < 90 with the plages. The specific heat 
was measured by an ac technique^ allowing us to mea- 
sure small samples (here ~ 100 nano- grammes). Heat 
was supplied to the sample at a frequency u> of the order 
of 70 Hz by a light emitting diode via an optical fiber. 
The induced temperature oscillations were measured by 
a 12/im-diam chromel-constantan thermocouple that was 
calibrated in situ. 



III. H-T PHASE DIAGRAM 



Fig. 1 displays the temperature dependence of the 
specific heat at various magnetic fields up to 7 T for 
-ff || a&- planes and up to 2 T for H perpendicular to those 
planes. In zero field, the specific heat discontinuity at 
T c , AC(T C ), represents only about 3 % of the total sig- 
nal, i.e. j-5 to 6 times less than in the best polycristallinc 
samplesu. However, our samples are extremely small and 
the addenda contribution are thus very large although 
impossible to estimate accurately. The data are there- 
fore given in arbitrary units and the temperature range 
is limited to T > 12 K. No significant increase of the 
transition width, AT ~ 2 K, is observed for (j,qH\\c < 
2 T and noH\\ab < 5 T. Despite a clear broadening 
above these chareristic fields, the specific heat anomaly 
remains well defined in the entire explored temperature 
range. These results imply that the broadening of the. 
transition width observed in polycristalline samplescm 
is due to randomly oriented anisotropic grains rather 
than to thermodynamic fluctuationsB. The curves are 
remarkably identical for both directions once the field is 
renormalized by the temperature dependent anisotropy 
coefficient (see inset of Fig. 1 for instance). It is also 
worth noticing that the amplitude of the specific heat 
jump AC decreases much too rapidly with H implying a 
rapid increase of the Sommerfeld coefficient 7 and there- 
fore of the density of states. An entropy conservation 
construction shows that 7 ~ 0.97jv (where 77V is related 
to the normal state density of states) for H ^0.5 H C 2 
in good agreement with .recent low temperature measur- 
ments by Bouquet et alU. 




10 15 20 25 30 35 40 

FIG. 1. Temperature dependence of the specific heat at 
different magnetic fields for H\\ab— planes (closed circles) and 
H± a b (open squares). In the inset : specific heat transition at 
[toH_i_ a b =1T and fj,oH\\ a t = r .fioH \Lab ~ 3.5 T showing that 
the same curve is obtained in both direction when the field is 
renormalized by F. 

The critical temperature has been determined using 
the classical entropy conservation construction. The cor- 
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responding H C 2 values as well as the one deduced from 
susceptibility (onset of the diamagnetic response) and 
transport measurments have been reported in Fig. 2. In 
the latter case, the critical field H C 2(T) or temperature 
T C (H) values have been defined at the onset of finite 
resistivity (see Fig. 3). 

Some of us have previously shown that the onset of finite 
resistivity coincides with the location of the specific heat 
anomaly for H±ab. Here we show that the same holds 
for parallel fields H\\ab. Note that the resistivity reaches 
the normal state value at a magnetic field Hr n higher 
than the as-deduced upper critical field (see Fig. 3). Be- 
low we argue that Hr n is due to the surface effects. 
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FIG. 2. H — T phase diagram of the MgB2 single crystal. 
Circles - H C 2 from transport measurments. Squares - H C 2 
from ac-susceptibility. Triangles - H C 2 from specific heat. In 
the insert : temperature dependence of the anisotropy. 

One can see that all three presented methods reveal 
the same results in the common range of applied fields 
(fioH <5 T). The upper critical field perpendicular to 
the basal plane has a typical temperature dependence 
for a type-II superconductor with a linear shape near 
the zero-field transition temperature and a saturation 
at the lowest temperatures with fioH C 2±ab — 3.5 T. On 
the other hand the parallel upper critical field has a dif- 
ferent strength and shape: close to T c , H C 2\\ a b reveals a 
positive curvature which changes to negative below 20 K 
and saturates to about 17 Tesla at the lowest tempera- 
tures. A direct consequence of the different forms of the 
temperature dependencies of these two critical fields is a 
temperature dependent anisotropy factor T. Then, T ^5 
is found at low temperatures but it is about 2 near T c 
(see inset of Fig. 2). 

The presented data extend the ones previously obtained 
by sonae of us for hqH < 7 T on crystals from the same 
batchli!]. It has been suggested that the positive cur- 
vature observed hece for _ff||a6 is a consequence of the 
two-gap structuretil. However, it is worth mentioning 
that a very similar behavior has also been observed 
in conventional superconductors including NbSe2c3, 
borocarbidesEJ and was believed to be a universal char- 



acteristic of layered compounds!^. This is even more 
"general" as it has,-also been observed in the isotropic 
(K,Ba)Bi03 systemli-1. The origin of this effect thus still 
has to be clarified. 
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FIG. 3. Magnetic field dependence of the MgB2 resistance 
at T = 5.4 K for various angles between the field and the 
ab— plane. The critical fieds H C 2 and H c $ are marked by the 
arrows for two principal field orientations. For these fields 
also the effect of current density is presented (dashed lines 
from right to left : j =17, 50, 170 A/cm 2 , solid lines j =500 
A/cm 2 ). 

Fig. 3 represents the resistive transitions from the 
superconducting to normal state taken at 5.4 K for dif- 
ferent angles between H and the a6-planes of the crystal. 
The measurments have been done at the current density 
500 A/cm 2 (solid lines). The current dependence of the 
transition for perpendicular fields (HJ-ab) - which is also 
highly non-ohmic - has been attributed to .surface con- 
ductivity in a previous paper by some of ust2l. We have 
found similar, albeit a bit weakerEj, effects for parallel 
fields (H\\ab) down to the lowest temperatures. In Fig. 
3 for both principal field orientations -ff||a6 and HJ^ab 
an effect of the measuring current is presented. It shows 
that the transitions are indeed current dependent. For 
current densities equal to 17, 50, 170 A/cm 2 transitions 
are smooth and for H±ab the onset of the finite resis- 
tivity is even shifted towards higher fields for decreasing 
current densities. For j > 500 A/cm 2 the onset of the 
transition is not changed anymore (measured up to 2000 
A/cm 2 ) and shows up as an abrupt jump at H = H C 2 
(for H\\ab already at 170 A/cm 2 ) followed by a smooth 
transition towards the normal state resistance Rn ■ The 
ratio between the fields at the onset of finite resistivity at 
500 A/cm 2 (Hr = o — H C 2) and the one at the end of the 
transition at the lowest current Hr n = H c3 is about 1.8 
for H-Lab and 1.4 for H\\ab. This ratio remains constant 
for all temperatures apart those very close to T c , where 
it is affected by the intrinsic width of the transition AT C . 
Our measurments can be compared to the ones obtained 
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in a pioneering work on surface effects in type-II super- 
conductors by Hempstead and Kimlij. They performed 
measurments on Nbo.sTao.5 and Pbo.83Ino.17 sheets. For 
fields parallel to the surface of the sheets (that is a 
favorable configuration for appearance of surface super- 
conductivity) they obtained results strikingly similar to 
ours. For sufficiently low current densities, the onset of 
a finite resistance is sensitive to the excitation current 
since the surface sheath could bear a current keeping the 
zero resistance even above H C 2- However, above some 
threshold current, this surface superconductivity is de- 
stroyed and the onset field does not change anymore; 
it has thus been identified with the upper critical field 
i/ c2 ll3. As obverved in our data, for H = iJ c2 the re- 
sistance of the sheets increased sharply to a fraction of 
Rn which was bigger for higher current and then, the 
resistance finally went smoothly to the normal state Rjy. 
H C 3 was best defined when measured at the lowest possi- 
ble current density, but the experimental H c3 /H C 2 ratio 
scattered between 1.6 and 1.96 and droped down to 1.14 
for copper coated Pb-In sheets. Indeed, the configuration 
of our transport experiment with H\\ab is very favorable 
for surface effects which must then be taken into ac- 
count. On the other hand, surface superconductivity 
effects should be negligible for HLab. However, in our 
experimental set-up the current and voltage electrodes 
overlapping the sample go over the vertical side planes 
which thus inevitably contribute to surface superconduc- 
tivity as well. 

H C 2 in MgB2 crystals has been measured by several au- 
thors in a limited temperature range. The H C 2± a b values 
of Angst et al.EZl obtained via torque magnetometry are 
in a good agreement with ours while -ff C 2|| a & arc slightly 
higher (implying a slightly higher anisotropy parame- 
ter). Similar H C 2r- values have also been obtained by 
Sologubenko et al.113 from the magnetic field and tem- 
perature dependencies of the thermal conductivity below 
6 Tesla. They also found a standard dependence for 
H C 2±ab(T), and a positive curvature for H C 2\\ a b{T) above 
30 K. Our data are in very good agreement with those 
obtained by the measurments of the magnetic moment 
in MgB2 single crystals by Zehetmayer et al£3. Finally, 
Bud'ko and CanfieldEll have arrived to similar conclu- 
sions by extracting the superconducting anisotropy from 
the magnetization measurments on randomly oriented 
powder samples. It was however of fundamental impor- 
tance to show that the "critical fields" obtained from 
all those measurments coincide with the thermodynamic 
transition deduced from specific heat measurments. In 
many novel superconductors fluctuations broadens the 
transitions from the superconducting to the normal state 
in magnetic field. One of the most striking phenomenon 
which has been observed in these systems is the existence 
of a melting line T m (H) above which the vortex lattice 
melts into a liquid of entangled lines. The presence of 
this liquid phase complicates the determination of the 
upper critical field H C 2- For instance, the onset of a finite 
resistance in those systems is indicating the melting tran- 



sition (at -fffl-o) while the upper critical field is located 
at a very end ofj-the resistive transition R{H C 2) = Rn 
(see for instancecil) . However, fluctuations are not ex- 
pected to play any significant role in MgB2 as confirmed 
by specific heat measurments: the width AC does not 
increase significantly with magnetic field and the upper 
critical field occurs near the onset of finite resistance at 
Hr—o. A non correct criterion for the H C 2 determination 
partly explains the discrepancy in the values deduced 
from transport measurmentsoEj. 
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FIG. 4. Angular dependence of the upper critical field 
(open circles) at 5.4 K together with the fits by the 
Lawrence-Doniach formula. Triangles are "critical fields" 
taken at the end of the resistivive transition Hr n showing 
the cusp- like behavior at small angles. In the inset: angular 
dependence of the critical current at fioH —2 T and T — 18 
K (Jo ~ 5xl0 3 A/cm" 3 ). 



Next, we discuss the angular dependence of the critical 
field. As shown in Fig. 3, the resistivity rapidly drops 
towards zero at hqH ~ 16 T for iJ||a& and progressively 
shifts towards lower fields as the angle increases. The 
corresponding H C 2(Q) can be well fitted by the the sim- 
ple ellipsoidal form ( g - 2 < e)sine ) 2 + ( H ^ e} cosQ ) 2 = 

1 of Lawrence and DoniachcJ for anisotropic three- 
dimensional superconductors with Y = 4.8 (see Fig. 4). 
A similar dependence has been observed at 26 K (not 
shown) but with V — 3.4. By triangles we have also 
plotted the characteristic fields Hr n . If surface conduc- 
tivity becomes important for iJ_La6 and ff||a& (Hr n is 
then close to ^3), this effect is diminished for other 
field orientations and the Hr n / H C 2 ratio drops down to 
~ 1.1. Surface conductivity effects can thus account for 
the cusp-like behavior of Hr n at small angles-idiich has 
been first attributed to H C 2 by other authorscijlHj. 



4 



IV. PEAK EFFECT AND ANGULAR 
DEPENDENCE OF THE CRITICAL CURRENT 

Another interesting phenomenon is the observation of 
a peak effect in the critical current. This effect shows up. 
in our transport measurments in the same way as int£l 
and is also visible in our susceptibility data in a small 
magnetic field range ~ 2-3 T for H perpendicular to the 
ab— plane. This increase in the critical current appears 
as a dip in the susceptibility in both temperature and 
magnetic field sweeps as alscv-jeported by Pissas et al. 
for the same field orientationc3. Surprizingly, the peak 
effect is only visible in a narrow magnetic field range. 
For lower magnetic fields, the peak is replaced by a rapid 
drop of the susceptibility and it has been suggested by 
Pissas et al.E3 that this drop could be associated with 
a first order transition in the vortex lattice in analogy 
with high T c oxides. However, in MgB 2 , the peak effects 
appears close the H c2 line and no sign for the existence 
of a vortex liquid could be obtained in our experiments. 
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FIG. 5. Magnetic field dependence of the critical current 
at HLab at T(K) = 8.5, 13, 17, 21.5 and 26 (from top to 
bottom) showing a peak effect at low temperatures (Jo ~ 5 
10 3 A/cm 3 ). In the inset : temperature dependence of the 
ac-suscetibility at [iqH —2 T for different angles between H 
and the ab— plane. 

The peak is here reflecting the fact that the shear en- 
ergy of the flux lattice drops towards zero more rapidly 
than the pinning energy in the vicinity of the normal 
state allowing to the vortex matter to accomodate more 
efficiently to the disorder close to the transition^. More- 
over, as shown in Fig. 5, the value of the peak current 
rapidly decreases for increasing temperature (i.e. for de- 
creasing magnetic field) but probably still exists at high 
temperature. Similarly, as shown in the inset of Fig. 



5, the dip in the susceptibility (that is the peak in the 
current density) also rapidly disappears as the magnetic 
fields is turned away from the c— axis. However, the peak 
effect is still visible in transport measurments for H\\ab 
again suggesting that the peak has shifted towards lower 
J values when the magnetic field was turned and thus dis- 
appeared from our experimental window. The J values 
presented in Fig. 5 have been deduced from the magnetic 
field dependence of the ac— susceptibility. Indeed, in the 
non linear regime, the ac response is related to the cur- 
rent density J through : x(Kc,u,T) = F(h ac / J(u/,T)) 
where F(x) depends on the pinning mechanism, the sam- 
ple geometry and the flux dynamics. F has been deter- 
mined experimentally following the procedure introduced 
by Pasquini et alEB and is in very good agreement-svith 
the numerical calculations performed by Brandtcj for 
bulk pinning in cylinders (with the appropriate thick- 
ness/radius ratio). Note that our critical current values 
are much smaller than the one previously obtained in 
polycrystalscj but, as shown below, J sharply increases 
when the field is applied parallel to the ab— planes. 
As shown in the inset of Fig. 5, at low temperature, the 
susceptibility (and thus the critical current) only weakly 
depends on 9 down to 6 ~ 10 degrees and the rapid 
drop of the susceptibility close to H c2 is progressively 
shifted towards higher temperature for decreasing 9 val- 
ues. Finally close to the ab— plane the susceptibility 
rapidly decreases showing that the pinning of the vor- 
tices is much more efficient for H \\ab. The corresponding 
J/ Jo values at 18 K are displayed in the inset of Fig. 4 
showing that J increases by a factor of about 4 for H 
parallel to the ab— planes. A similar behavior has been 
observed in the critical current deduced frjQm transport 
measurments in c— axis oriented thin filmsEH and Eltscv 
et al. recently reported on a rapid drop of the dissipa- 
tion when the magnetic field-is applied parallel to the 
ab— planes in single crystalsu. A similar cusp in the 
critical current has been observed in high T c cuprates in 
which the vortices are strongly pinning by the weakly su- 
perconducting layers between the Cu02 planes when the 
field is applied parallel to those planes (so called intrinsic 
pinningO). However, the origin of the strong pinning of 
the ab— planes in diborides still has to be clarified given 
the rather small anisotropy of this material. 



V. CONCLUSION 

In summary, the upper critical field has been deduced 
from specific heat, high magnetic field transport and 
ac-susceptibility measurments for H parallel and per- 
pendicular to the ab— planes. H C 2 perpendicular to the 
planes reveals a conventional temperature dependence 
with (J<oH C 2±ab(0) — 3.5 Tesla but the parallel critical 
field with HoFt C 2\\ab{fy — 17 Tesla has a positive cur- 
vature at temperatures above 20 K. Consequently, the 
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anisotropy factor T = H C 2\\ a b/ Hci^ab is temperature de- 
pendent decreasing from ~ 5 at low temperature towards 
~ 2 close to T c . The critical current deduced from the 
susceptibility measurments present a sharp peak effect 
for noH-Lab ~ 2 — 3 T. At low temperatures, J only de- 
pends weakly on the angle between H and the ab— plane 
for 9 > 10 degrees but rapidly increases for H\\ab. 
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